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Abstract 

The Kumaraswamy distribution has been proposed as an alterna- 
tive to the beta distribution with more benign algebraic properties. 
They have the same two parameters, the same support and qualita- 
tively similar shape for any parameter values. There is a generic re- 
lationship between the distributions established by a simple transfor- 
mation between arbitrary Kumaraswamy-distributed random variables 
and certain beta-distributed random variables. Here, a different rela- 
tionship is established by means of the logistic map, the paradigmatic 
example of a discrete non-linear dynamical system. 
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1 Introduction 
1.1 The Logistic Map 

The logistic map is a parametric discrete-time non-linear dynamical system. 
It is widely studied because of its complex transition from o rderly to ch aotic 



behavior in spite of an extremely simple defining equation [Mavl . Il976l | . 



Definition 1 (Logistic map). For a real parameter r > 0, the logistic map 
is the function 

f r (x) = rx(l — x) (1) 
restricted to the closed real interval [0, 1]. This is a total function for r < 4. 

The objects of interest are the trajectories that arise from the iteration 
of f r , i.e., sequences with the recurrence relation = f r (x n ). With in- 

creasing r, they show all kinds of behavior from convergence in a singular 
attr actor through bifurcating periodic solutions to deterministic chaos. The 
case r = 4 is of special interest: It is known to be chaotic and ergodic on the 
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interval [0, 1]. The density of the attractor is given by the beta distribution 
with parameters a = f3 = 1/2 and hence, by ergodicity, the long-term dis- 
tribution of states under almost all initial conditions (with the exception o f 
unstable periodic solutions) converges to that distribution [Jakobsonl . Il981 |. 



1.2 The Kumaraswamy Distribution 

Kumar aswamvl 



1980l | introduced his now eponymous distribution, originally 
called double-bounded distribution, as an alternative to the beta distribution. 
They have the same real parameters a, f3 > 0, the same support and similar 
shapes, but the Kumaraswamy distribution function, unlike the beta distri- 
bution function, has a closed algebraic form. Is has been found both more 



accurately fitting hydrological dat a in s i mulat ions [Kumar aswamvl . Il980l | and 



computationally more tractable Jonesl . 12009], 



The similarity between the two classes can be formalized. It is known 
and easy to see (from equations (|3|) and (|6|) below) that, if random variable 
X is Kumaraswamy-distributed with parameters a and /?, then X a is beta- 
distributed with a = 1 and the same j3. Here, a new formal relationship 
between certain beta and Kumaraswamy distributions is established. 



2 Computing State Distributions of the Logistic Map 

Given a random variable X with known distribution on [0,1], the propagated 
distribution of f r (X) can be computed. In the following, the term continuous 
probability distribution (cdf) refers to a cumulative distribution function with 
support [0, 1], that is, a continuous weakly monotonic function F : [0, 1] — > 
[0, 1] with -F(O) = and F(l) = 1. A random variable X is distributed 
according to F, written X ~ F, if and only if P{X < y) = P(X < y) = F(y). 

Definition 2 (Propagation). Let F be a cdf. Then f r (F) is another cdf, 
namely 

f r (F)(y) = F(l-q r (y))+l-F(± + q r (y)), (2a) 

where 

= (2b) 
[0 ify>|. 

Theorem 1. If X ~ F then f r (X) ~ f r (F). 
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Proof. 



P(f r (X) <y) = P{rx(l -x)<y) 

= P{\x-\\>q r {y)) 



P(x< \-q r {y)yx> \ + q r {y)) 
P{x<\- q r {y)) + P{x > \ + q r (yj) 
P{x<\- q r (y)) + 1 - P(x < \ + q r (y)) 
F{\-q r {y))+l-F(\+q r {y)) 

Jr{F)(y) 



one finds that 



h(A){y) = 2A{\ - q A (y)) = 2A ( ± = A{y) . 



□ 



The beta distribution with real parameters a, (3 > has the cdf 

B{a,(3){y) = £ ^(1 - tf^dt / £ ^(1 - tf~ l dt . (3) 

For a = /3 = 1/2 the beta distribution is the arcsine distribution. 

5(1/2, l/2)(y) = A(y) = 2/vr arcsin (4) 

Theorem 2. The arcsine distribution is a fixed point of f^. 
Proof. Using the symmetry of (HJ) and the trigonometric identity 



_ . 1 — y/1 — y 
arcsin yjy = 2 arcsin y , (5) 



□ 



The uniform distribution has the cdf U(y) = y. The Kumaraswamy 
distribution with real parameters a, (3 > has the cdf 

K(a,(3)(y) = l-(l-y a f. (6) 

Theorem 3. Let X be distributed uniformly. Then f/^(X) = f&if^X)) is 
Kumaraswamy- distributed with a = [3 = 1/2. 

U\u) = K{\,r) (7) 
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Proof. For r = 4, (I2bp simplifies to (74(2/) = a/1 — y/2 and hence (j2a| to 



Then after one step one has 



such that fi{U) = K(l, 1/2), and after two steps one has 
h 2 (U)(y) = h{h(U))(y) 



The term s looks rather convoluted, but its square simplifies to s 2 = 1 — -^/y, 
so with s > conclude that s = ^1 — ^fy. Hence from ([6j follows (|7J): 

/ 4 2 (t/)( y ) = 1 - JlZjjj =!-(!- y 1 / 2 ) 1 ^ 



and thus fi 2 (U) = K (1/2, 1/2). □ 



3 Conclusion 

We have examined the probability distribution of successive states of the 
chaotic dynamical system described by the logistic map with parameter value 
r = 4. It is known that this system is ergodic and that its long-term behavior 
is described by the beta(l/2, 1/2) distribution. Hence the distribution of 
states will converge almost surely. In order to investigate individual elements 
of that convergent sequence, we have given a closed formula for propagating 
probability distributions over the iteration of f r (Theorem [T|, and validated 
this approach by demonstrating that the beta(l/2, 1/2) distribution is indeed 
a fixed point solution (Theorem [2]) . We have proceeded to demonstrate that, 
starting from the uniform distribution (the maximum-entropy distribution on 
the state space), after two steps the Kumaraswamy(l/2, 1/2) distribution is 
reached (Theorem |3J). Figure Q]shows the cdfs of the initial state and the first 
four iterates plotted over the uniform, Kumaraswamy and beta distributions. 

The above findings qualify as purely theoretical results in the sense that 
they provide original and provable insights about the relationship of math- 
ematical objects of theoretical interest. No suggestion can be given yet as 
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to the application to real- world statistical problems. Apart from the investi- 
gation of potential practical implications, the presentation of higher iterates 
f^{U), ... is an open problem. The initial sequence of distributions fi n {U) 
looks pretty regular: 

U = B{l,l) = K{l,l\ —> K(l,l) —> K{1\), 

n=0 n =l n=2 

but there is no obvious continuation of the pattern for n > 3. 
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Cumulative Distributions 




Figure 1: State cdfs D n = f^ n {U) of an initially uniformly distributed sys- 
tem running the logistic map with r = 4, after n = 0, . . . , 4 steps; plotted 
over the uniform (C7), Kumaraswamy {K) and beta [B) distributions, with 
parameters a = f3 = 1/2. D$ coincides with U. D2 coincides with K. For 
n > 3, D n converges towards B. 
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